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j^ ' Abstract 

r^ ■ In many problems of quantum chaos the calculation of sums of 

products of periodic orbit contributions is required. A general method 
of computation of these sums is proposed for generic integrable models 

^ ' where the summation over periodic orbits is reduced to the summa- 

tion over integer vectors uniquely associated with periodic orbits. It 
is demonstrated that in multiple sums over such integer vectors there 
exist hidden saddle points which permit explicit evaluation of these 
sums. Saddle point manifolds consist of periodic orbits vectors which 
are almost mutually parallel. Different problems has been treated by 
this saddle point method, e.g. Berry's bootstrap relations, mean val- 
ues of Green function products etc. In particular, it is obtained that 
suitably defined 2-point correlation form-factor for periodic orbit ac- 
tions in generic integrable models is proportional to quantum density 
of states and has peaks at quantum eigenenergies. 
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1 Introduction 

The trace formulas are the main instrument in the investigation of relations 
between quantum and classical properties of a given system (see e.g. [|l|]). 
These formulas connect the quantum density of states, d{E), in the semi- 
classical limit with a sum over classical periodic orbits (po) ^-^ 

d{E) = d{E) +J2A,expit^^), (1) 

po '^ 

where d{E) is the mean density, Sp{E) is the classical action along a periodic 
orbit, and Ap is a pre-factor build from classical quantities. 

There is two main applications of these formulas. First, one tries to use 
them to compute smoothed density of states (or even approximate positions 
of energy levels 0). Second, one uses trace formulas to understand statistical 
properties of energy levels [^ . In the latter approach one starts with a formal 
expression of n-point correlation function, i?„(ei, e2, . . . , e„), as a mean value 
of product of n density of states, d{E), at different points 

i?„(ei, €2, . . . , en) =< d{E + e,)d{E + ea) ■ ■ ■ d{E + e„) >, (2) 

where the brackets denote an average over an interval of energy AE which is 
classically small but includes a large number of quantum levels. Substituting 
in this formula the semi-classical expression ([I|) for the density of states one 
obtains a semi-classical approximation for correlation functions. 

The main difficulty of such method is the calculation of mean values of 
products of contributions of different periodic orbits 

< E exp(4^e,5p,(i?))>, (3) 

Pl,P2,--- j 

where Cj = ±1. When all Cj are of the same sign this average equals zero 
but if signs of Cj are different the calculations are far from being clear. 

The simplest approximation (called the diagonal approximation) consists 
of taking into account only terms with exactly the same action 

<exp4(S„-S„J>={j;;;|;| . (4) 



Berry [^ showed that for integrable systems this approximation gives the 
correct 2-point correlation form-factor but for chaotic systems the vahdity of 
diagonal approximation is limited to short - time behavior of the form-factor. 
For more complicated quantities (even for integrable systems) this approach 
is insufficient and more refined methods are needed 01-[1I- 

The purpose of this paper is to develop a method which permits explicit 
computations of the expressions (|^) for generic integrable systems. We shall 
show that in these expressions there exist hidden saddle points which give 
dominant contribution in the semi-classical limit. 

The usual trace formulas can be applied only to quantum mechanical 
quantities which (like in Eq. (|l])) are expressed through one Green function. 
The method discussed in this paper permits to construct a new type of trace 
formulas which express a product of such quantities for integrable systems 
through a sum over classical periodic orbits. 

The plan of the paper is the following. In Section 2 for completeness the 
derivation of the trace formula for a rectangular billiard is presented. This is 
a simple prototype of integrable models where all formulas are transparent 
without unnecessary complications. In Section 3 on the example of Berry's 
bootstrap we demonstrate the existence and importance of quantities which 
include sums like in Eq. (^) and whose semi-classical calculation is the main 
topic of the rest of the paper. The saddle point method for oscillatory and 
smooth terms is discussed in Sections 4 and 5 respectively. In is shown that 
the dominant saddle point configurations are build from periodic orbits whose 
integer vectors are almost mutually parallel. In Section 6 the calculation of 
mean values of products of different powers of retarded and advanced Green 
functions is presented. In this Section we also discuss the calculation of next- 
to-leading terms by the saddle point method. The correlation functions for 
classical actions are derived in Section 7. An interesting consequence of the 
discussed formalism is the fact that (properly defined) 2-point correlation 
form-factor for actions of classical periodic orbits is proportional to quantum 
density of states. The generalization of the proposed method for generic 
integrable models is performed in Section 8. 



2 Trace formula for rectangular billiard 

Let us consider a plane rectangular billiard with sides a and h with (for 
simplicity) periodic boundary conditions. The functions 

, , 1 .2-111 2711 , . , 

'^m,n[^, y) = —^ exp( nx + ——my), (5) 

yah a b 

with integers m, n = 0, ±1, ±2, . . . are periodic solutions of the (Schrodinger) 
equation 

(E„„ + A)^^,„(x,i/)=0 (6) 

with energy 

Emu = i-fn' + (^)W. (7) 

a 

The quantum density of these states 

oo 

d{E)= Y. SiE-E„J (8) 



m,n=— oo 

may be rewritten by the Poisson summation formula in the following way 

,27r,^ o .27r, 

M,N=-oo ' 



d{E)= g /e2-(^^™+^")(5(E-(— )V-(^)W)dncim. (9) 



The substitution n = ar cos 6/{27i) and m = br sin ^/(27r) after simple algebra 
gives 

c^(E) = - ^ Jo(A:^Miv), (10) 

'*'"" M,Ar=-oo 

where E = k"^. 



Lmn = ^{NaY + (M6)2, (11) 

and 

Jo(a;) = — /'"e'"=°^^rf^ (12) 

27r JO 

is the Bessel function of zero order. It is evident that Lmn is the length of 
a classical periodic orbit or, more strictly, the length of a classical trajectory 
on a resonant (= periodic) torus. In integrable models almost all periodic 



orbits belong to resonant tori and we shall use these two notions on the equal 
footing. 

As most of the terms in both parts of Eq. (|TUp are 4 times degenerated it 
is convenient to define the density of non-degenerate states by dividing both 
parts of this equation by 4 

d{E)=d{E) + (r"'{E). (13) 

Here 

m ^ A. (14) 

where A = ah is the area of the rectangle and 

d°^\E) = ^Y.^,J,{kL,), (15) 

Lip 

where the sum is taken over all periodic orbits defined by two integers M 
and A^ whose lengths Lp are given by Eq. (pi]), e^ = 1 if both M and N are 
nonzero and e^ = 1/2 otherwise (the term with M = N = {] gives d{E)). 
In semi-classical limit k —>■ oo and 

Eqs. ([T3|)-(p!6[) define the trace formulas for a rectangular billiard (with pe- 
riodic boundary conditions). Trace formula for general integrable systems 
have similar form and will be discussed in Section 8. 



3 Berry's bootstrap 

The above trace formula for rectangular billiard is the simplest example 
of trace formulas for integrable systems. We know explicitly everything. 
The periodic orbit lengths are given by simple expression (|1T]) and all semi- 
classical corrections to Eq. (|l^) (coming from asymptotic expansion of the 
Bessel function) are easy to take into account. 

But even in this case quantities like those in Eq. (|^) is difficult to compute. 
The necessity of such calculation is clearly seen e.g. from Berry's remark [H 



that for any quantum system with non-degenerated discrete spectrum the 
square of the density of states should be proportional to the density itself. 

This fact can be demonstrated as follows. The density of states at real 
values of energy is usually defined as the limit e — » of the sum over all 
energy eigenvalues e„ 

"^'^^^ ^2m^^E-en-te~E-en + te^^n^iE- e^f + e^ ' ^^'^^ 

In the semi-classical trace formula a finite value of e corresponds to a regu- 
larization of divergent sums by adding (e.g. in Eq. (0)) a small imaginary 
part to the energy E —>■ E + ie. 

As it is evident from Eq. (0) at finite (but small) values of e the function 
d^{E) has peaks near each energy eigenvalues, the pear width being propor- 
tional to e. Now let us compute d'^{E) assuming that all eigenvalues are 
non-degenerated. It is clear that in the limit e — > 

and cross terms with different En will give negligible contributions at small 
e. 



Because 



one gets that 



and therefore 



dx vr ^^g^ 



oo 



(x2 + e2)2 2e3' 



e^ 71 



limrr— ^^= o-^l^), (20) 



\im27red^^{E) = d{E). (21) 

Similarly by computing higher powers of d^{E) one obtains that 

These bootstrap relations reflect the analytical structure of the density of 
states, namely that it should have (5-function singularities at real values of 



energy with unit residues. Therefore these type of relations is very general 
and our first purpose is to find a method which will permit to prove these 
relations for integrable systems starting from the semi-classical trace formula 

4 Saddle points for oscillatory terms 

Let us begin with the case n = 2. The density of states may be written as a 
sum of two terms de{E) = d+^^E) + d^^{E) where ± corresponds to the sign 
of the exponent in two terms in Eq. (p!6D 

rl , r/?l - V ^^ ikLj,-i-K/i-eLj,/2k /^o^ 



and d_^{E) = d\^{E). Of course, the contribution of d{E) to the left- 
hand side of bootstrap relations is negligible at small e. Therefore d1{E) = 
d\^{E) + 2d+,{E)d_,{E) + d'i^{E). The periodic orbit expansion of d'^^^{E) 
contains only terms with the same sign of the actions in the exponent and for 
any fixed total action it includes only a finite number of terms, and in the 
limit e — > their contributions are negligible. On the contrary the periodic 
orbit expansion of d^^{E)d-^{E) for any finite total action contains infinite 
number of terms and it is this term which can give contribution in the limit 
of small e. Therefore the following relation should be true 



lim27re — :r T^ 



^ -^ ^Pl'^P2 ^ik(h-l2)-(li+l2)e/2k 

27r2 



Pl:P2 




e 



d{E) + — V( /^ eiH,^»^/4 ^ ^_^y (24) 

2vr p JSnklp 

Here in the left-hand side li^2 are the lengths of periodic orbits pi and p2 and 
the sum is taken over all such orbits. 

The smooth part (=mean value) of this expression is easy to compute. 
It was shown in that diagonal approximation is applicable and the mean 
value of the level density is expressed through the sum over periodic orbits 
as follows 

42 p-d/k 

lim27re^-y--- = d(E). (25) 



The density of periodic orbits can easily be calculated from Eq. (0). If N{1) 
is the number of periodic orbits with the lengths less than / then for large / 

Nil) = ^. (26) 

where A = ab is the area of the rectangle and the factor 1/4 comes because 
we consider only orbits with positive integers (M, A^) . 

Changing the summation over periodic orbits to the integration over pe- 
riodic orbit length with the above density one concludes that the left hand 
side of Eq. ( P5[ ) in the limit of small e equals 

2ne — - / -—Idl = . (27) 

2n^ Jo IGAkl Wn ^ ' 

But this is exactly equal to d{E) as it was predicted by the bootstrap equation 
Though in this case the computation of the smooth part is simple the 



bootstrap formula (|2T| ) should be valid for oscillatory terms as well for which 
the diagonal approximation cannot be applied. It is the calculation of os- 
cillatory terms (proportional to exp{iklp)) which is the main subject of this 
Section. 

The left-hand side of ( PT| ) contains a double sum over pairs of periodic 
orbits. Let us try to find saddle points in this double sum. 

A periodic orbit for a rectangular billiard (and for any integrable models 
in 2 dimensions (see Section 8)) is defined by two integers M and N (or integer 
vector (M, A^)) and its length equals the modulus of vector with components 



Ma and Nb, L{M,N) = ^J{May + {Nby. The length of a periodic orbit 
defined by integers M + 6M and A^ + SN has the following expansion up to 
the second order in 6M and 6N 

L{M + 6M,N + 6N) = ^a^{M + dMf + b'^{N + 5Nf (28) 

r,.. .r^ a^M6M + b^N6N A^ ,^^„^ ,....^2 



The double sum in Eq. ( pi]) has the form Y.pi,p2 ^W ^^{li — I2). Each peri- 
odic orbit in this sum is defined by two integers and the sum is really the 
sum over 4 integers. Let us assume that saddle point manifolds consist of 



periodic orbits with integers (Mi + 5Mi,Ni + 5Ni) for the first sum and 
(M2 + (5M2, N2 + 5N2) for the second one with unknown integers (Mj, iVj), 
(5Mj,(5A^j). The necessary condition for the existence of a saddle point in 
the double sum is the cancellation of linear terms in the exponent. From 
Eq. (p8[) one concludes that the saddle point condition has the form 

a'^MiSMi + b'^Ni5Ni a^M26M2 + b'^N25N2 

L{Mi,Ni) L{M2,N2) ~ ■ ^ ^ 

The important restriction is that we are looking at the integer solutions of 
this equation. We shall also assume that the ratio a^/b"^ is a 'good' irrational 
number (which is a necessary condition that the spectral statistics of such 
billiard will be close to the Poisson statistics [10|). In this case the only 



possibility to find non-trivial solutions of the above saddle point equation is 
to require that the lengths L{Mi,Ni) and L{M2,N2) are commensurable 

L{M2,N2) ra' ^ ^ 

with certain integers ri and r2. This condition means that saddle point values 
of integers Mi and Ni are the following 

Ml = rim, M2 = r2m, Ni = rin, N2 = r2n, (31) 

where the pair of integers m and n has no common factors (i.e. they are 
co-prime integers). 

Now the saddle point condition (^9]) takes the form 

a^mSMi + b'^n6Ni - a^m6M2 - b'^n6N2 = 0. (32) 

As we have assumed that the ratio a^/6^ is an irrational number the only way 
to fulfill this equation is to cancel terms in front of a^ and 6^ separatively 

6M1 - 6M2 = 0, 6N1 - 6N2 = 0. (33) 

This argumentation shows that in the double sum ( pi] ) there exists saddle 
point manifold consisted from periodic orbits defined by the following pairs 
of integers 

Ml = nm + 5M, A^i = Tin + 5N, M2 = r2m + 5M, N2 = ran + SN, (34) 



for arbitrary co-prime integers m and n and arbitrary 6M and 6N. 

Note that the knowledge of pairs (Mi, A^i) and (M2, A^2) uniquely defines 
the pair {m, n) and the difference ri — r2 provided that it is non-zero. Namely, 
Ti — r2 is the greatest common factor of the pair (Mi — M2, Ni — N2) and 
{m, n) is the ratio of the division of (Mi — M2, A''i — A^2) on ri — r2. 

To find the value of the double sum (pT]) in the saddle point approxima- 
tion it is necessary to perform the summation over all saddle point manifold 
defined in Eq. (|3^ in the Gaussian approximation (i.e. expanding all actions 
up to the quadratic terms in 6M and SN as in Eq. (|28|) ). 

Let us denote the left-hand part of Eq. ( pi]) by D{E). Using ( p8D and 
lif) one obtains 



exp{ik{ri - r2)k - (n +^2)77^ + ^k—^{- — )^^) + c.c, (35) 



8^'^ io7^,r2 smN loVn^ 

2fc^ 2/3^1 r2' 

where t = nSM — mSN, and /q = va?rn?~+lPn^ . 

Note that the quadratic form in the exponent is not positive definite. But 
it is easy to check that if the quantity t has two identical values 

nSMi - m6Ni = n5M2 - m6N2, 

and integers m and n have no common factor (as it was assumed above) 
then 6M2 = 6M1 + Im and 5N2 = SNi + In for a certain integer /. But 
from Eq. ( |5^ it follows that this analog of zero modes corresponds to a 
change of the repetition number: r2 = r2 + /. Therefore the restriction of 
the summation to orbits with fixed repetition numbers is equivalent to the 
summation over all values of t only once. 
It means 

A^e °° 1 el 

X t <'M'k -^)jL~J''' t') + C.C. (36) 

t=-oo ^h^ir2 

Let fi = r2 + r. Changing the summation over t to the integration one finds 



D{E) = _ii_ V V^L^exp(^fcr/o-r^-^-) V e-'^^'^" + c.c. (37) 

10 



At small e the last sum tends to k/{elo) and 

A °° 1 e vr 

«(^) = ^ E g ^S^^ -M'krl, - r- - ,-) + C.C. (38) 

which is exactly the right hand side of Eq. ([2l|). 

Note that the changing the summation over t to the integration in the sum 
X]i^_oo^^P(^^^'^^) is equivalent to ignoring terms of the form e"*'^"* 1^ j ^Jx 
with integer m 7^ 0. This is correct when x ^ and only in a weak sence 
(i.e. after the multiplication of both sides of the equality 

^ giTTXt^ _ ^ _ ^ ^ ^-^.va^|. (3g) 

t=-oo V ^ V ^ m=-oo, m^O 

by a suitable chosen test function). It is in such weak sence that one should 
treat equalities similar to Eq. (^). This is not a restriction of our method as 
all semi-classical trace formulas have mathematical meaning only in a weak 
sence. 

The above calculations clearly demonstrate that the proposed saddle 
point method is sufficient to obtain the bootstrap condition pT| ) for n = 1. 
Below we show that exactly the same considerations permit to verify the 
bootstrap conditions (p^ for all n. 

Because dJyE) = d^^{E) + d_^{E), d'^{E) is given by the sum 

d:{E) = Y: C:^d'i\iE)dZiE), (40) 

where I'i > and z/i + z/2 = n. Substituting instead oid±^{E) its semi-classical 
expression ([l6|) one gets 

A ik{L-L')-e{L+L')/2k-in{ui-U2)/4: 

dZ{E)d''\{E) = ( Y^+-^ V , , (41) 

and L = li + . . . + li,-^, L' = l[ + . . . + l^^. The summation here is performed 
over all periodic orbits with lengths Ij, j = 1, . . . , z/i corresponding to terms 
with positive exponent and over orbits with lengths I'f., k = 1, ... ,1^2 from 
terms with negative exponent. 

Let each 'positive' periodic orbit is defined by a pair of integers {Mj, Nj) 
and respectively the 'negative' orbit is associated with a pair {Ml.,N'j^). The 

11 



same arguments as above prove that the saddle point manifold (i.e. the set 
of integer vectors for which linear terms in the difference L — L' cancel) has 
the following form 

Mj = rjm + 6Mj, Nj = VjU + dNj, M^ = r'f^m + 6Ml^, N'j^ = r[n + 6Nl^, (42) 



with 2 restrictions 

j=l k=l j=l k=l 

It is convenient to rewrite these conditions in terms of integer vectors 



(43) 



associated with each periodic trajectory. Now Eqs. 



(44) 
(H3|) take the form 



vi 



V2 



3=1 k=l 



rn. 



where 



and 



vi 



VI 



j=i k=\ 



(45) 



(46) 



(47) 



n = (m,n). 

Because we assume that m and n are co-prime integers they are defined 
uniquely from the knowledge of Nj and A^^ and each term in the multiple 
sum in Eq. ( ^If) can be uniquely attributed to the new summation of the 
form of Eqs. (|2|)-(|^). 

The total contribution of these saddle points to the n-fold sum in Eq. (p]) 
in the Gaussian approximation is 






where 



P 



A 



27rv/8^fH^' 



(48) 
(49) 
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and 

(50) 
Here 

Q = k^^, (51) 

R = ri + ... + r^^,R' = r[ + ... + rl^,tj = nSMj-mSNj, t'^ = nSM'^-mSNi 
and the (5-function is a consequence of the saddle point conditions ( ^31) in 
variables t and t'. 

Changing the summation over t and t' into the integration and using the 
usual representation of the ^-function 

5{x) = — fdae'"'', (52) 

one can easily perform all integrals and one obtains 

\l \r\ Q ' 

where r = R — R' (note that for oscillating terms r 7^ 0). 
This results leads to the following expression for dJ^{E) 

d:{E)=Y.P\h,T-' E N{R,,R,f- = , (54) 

/o,r V ^ R^-R2=r \J \r\ 

where 

n 

N{R,, R2) = E C:'n{R,, R2; z/i), (55) 

and n(i?i, R2; z/i) is the number of terms with fixed number of positive, -Ri, 
and negative, R2, repetitions. 

To find this number let us first compute the number of representation of 
a given number M into a sum of k integers: M = rii + . . . + n^ and all Ui > 1. 
It is easy to see that this number is the coefficient of x'^ in the expansion 
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x^ /{I — xY^ therefore there exists C\j\ ways of representing M as a sum of 
k non-zero integers and 

n— 1 

N{R,, R,) = Y: C:^Cgz\C-^;::\-\ (56) 

1/1=1 

Because we are interested in the region Ri,R2 -^ oo with fixed difference 
R1 — R2 = r we can use the following asymptotics of the binomial coefficients 



jj^i-i 






Finally 

iV(i?i, R,) = Y: C7C:\','—^ (57) 

The remaining sum 

n— 1 n— 1 

E/^ui/^ui — 1 \ ^ /^ui/^n—ui — 1 

^n ^n~2 — 2^ ^n ^n-2 
ui=l v\=l 

is the coefficient of x""^ in the expansion (1 + a;)"(l + a;)"~^ therefore 

n-l 



The total contribution is 



^ gifcZo»--J7rSgn(r)/4 oo J^'^"^ 



/o,r ^ \/ r R2=l '^"' ^''• 



When e ^ the last sum tends to {k/el^Y ^.nd because 

one obtains the following relation 

hm ^--l^d:{E) = -J: ^== + c.c. = d-{E). (61) 

14 



As 

4"-^ ^ {2n-2)\\ 

C--\ (2n-3)!! ^ ^ 

this results coincides with oscillating part of the general bootstrap condition 
(H) for all n. 

These calculations demonstrate that our saddle point method reproduces 
correctly the oscillating part of bootstrap condition but it remains the ques- 
tion how to apply this method for the smooth part of the bootstrap condition 
with n> 1. 



5 Saddle points for smooth terms 

Immediate difficulty of generalization of the above method to the computa- 
tion of the smooth part of the bootstrap condition (^) in the fact that now 
one cannot attribute uniquely values of (m, n) to each term in the multiple 



sum over periodic orbits (^I]) . Even for n = 2 the calculation of the smooth 
term is done by the diagonal approximation and requires the knowledge of the 
number of periodic orbits in contrast to the saddle point method discussed 
in the above Section. 

Each periodic orbit of rectangular billiard is defined by an unique vectorQ 

N = {M,N), (63) 

with positive integers M, iV and the length of this periodic orbit is given by 
the following expression 



L{N) = ^{aMf + {hNf. (64) 

In the polar coordinates with < < 7r/2 

M = Rcos(j), N = Rsincf), (65) 

the local density of periodic orbits p{l, (p) may be computed from the relation 

p{l, (t))dld<t) ^ / dMdN5{L{N) - l)dl, (66) 

^For the purpose of this Section it will be natural to consider vectors with component 
(Ma, Nb) but to be consistent with more general case discussed in Section this definion 
is more convenient. 

15 



which le 


;ads to 


n(] M — 


and 






-^0(0) = \l ci^ cos^ + 6^ sin^ 0, 



(67) 
(68) 



^/2 d<p ^ _T^ 

Ll{<P) 2ab' ^^ 



the integrated density of periodic orbit length 

r/2 ^, ,_ , nl 



p{i)=[ /^(^'<^)^<^=^' (70) 



in agreement with Eq. (pG]). 

To apply the saddle point method it is necessary to construct a config- 
urations of vectors ( |63D (saddle point manifold) such that any variations of 
them will decrease (or increase) the total length (i.e. the expansion of the 
total length will not contain terms linear on variations). From geometrical 
considerations it is clear that the only possibility for such a configuration is 
the case when all vectors are parallel. To build the saddle point manifold we 
shall proceed as follows. 

Let us consider n vectors A^j i = 1, . . . , n in polar coordinates 

Ni = Ri{cos(f)i,sm(j)i). (71) 

The condition that all these vectors are almost parallel is equivalent to the 
statement that all polar angles are close to each other 

0i = + 50i, (72) 

and 6(j)i <^ (p. Under these conditions L{Ni) can be calculated from Eq. (^) 
up to the second order of 6(f)i 

L{N,) = R^Lo + ^^^^1^}^^R,6<P, + ^R^{6<p.f. (73) 

ZLq ZLq 

Let us now compute the difference of actions 

n n 

Sr. = k{J2L{N,)-L{J2Ni)). (74) 

1=1 i=l 
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Note that this operation is possible because the sum of integer vectors (|6^) 
with positive components is also an integer vector with positive components 
and therefore it defines a periodic orbit. Direct application of Eq. (|T5D gives 



i=l Ki + ... + Kn i^^ 

= -^—^-—^T.RiRM-<P^?, (75) 



where 






As in the previous Section 



d:iE)= y: c:^dziE)dziE), 

ui+i'2=n 

and under the same notation as above 



exp(zfc(L-L')-^^^^^-.^(z/i-z.2)), (77) 

with L = li + ... + l^^ and L' = l[ + . . . + l'^^. 

Each periodic orbit is defined by an integer vector Ni and according to 
the above discussion the saddle point manifold for the smooth part of this 
expression is constructed from vectors Ni and N', which are almost parallel 
each to other and obey the following relation (because it should be the smooth 
part) 

Y.N, = T.Ni (78) 

j=i k=i 

One may easily check that the regrouping terms in Eq. (^) according to the 
sum of vectors Nj and that of N^ avoids the double counting and permits 
the correct arrangement of different contributions. 

In polar coordinates the saddle point condition ( [78| ) means that 

i?i + ... + i?,, = i?; + ... + K^, (79) 
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and 

Ri5(t)i + ... + RuMur = R'M[ + ... + KM'u2- (80) 

From general considerations (and from Eq. (fTSD) it follows that the difference 
of lengths {L — V) depends only on differences of angles between different 
vectors and therefore one of these angles (say 0i = 0) can be chosen arbitrary 
and under the condition of smallness of these differences the sum over periodic 
orbits in Eq. ( |77D can be substitute by the integration over components of 
Ni and Eq. ( ffTl) can be transformed as follows 

j-tt/2 c '^i '^a t^i i"i 

d'2,iE)dZiE) = / ^p^^^^'^^) / n dRj n dK n d<p, n d<p'>. x 

■^° '' j=l fe=l j=2 fc=l 



6{R - R')6{R<f) - R'<p')^{Ri ■ ■■Ru,){R[ ■ ■ ■ K,) x (81) 

exp(^Q(|: P,(0,)' - E KW) - ^^ \f^' - 4(^1 - ^^2))' 

j=2 k=l 



where R = Ri + . . . + R^„ R' = R[ + . . . + R'^^, Rcj) = i?202 + • • • + Pz.i0^i, 

(82) 



R'(f)' = R[(j)[ + ... + Rl^(f)l^, and 



1/2' 
1/2 "^ 1^2' ^^^ 

2n^8TrkLoi(l)) 



Note the absence of 0i in these expressions. 

Using the standard representation of the 5-function 

6{R(I) - R'cj)') = ^ r e*"(^'^-^'<^'), (83) 

ZTT J-oo 

one can easily perform the angular integration over all independent angles 

(j)j and 0'fc 

J7r(!/i-!^2)/4 

< ^e*(^-^ )5(P0 - R!<P') >= rZLv-i+-2)/2-i ^ 



(84) 
Combining this result together with Eq. ( |77D one obtains 

#p(-+-)(-, 
/o Q 

poo ^1 ^2 

/ n dRj n ^^fc-^i^ - p>-(«+^')^o/2^ (85) 
•^° j=i fc=i 



< drAE)drAE) >= r rf0p(-i+-2)rZL^(-i+-2)/2-i 

JO 



The integrals over R and R' may be computed either by using the repre- 
sentation (|83|) of the ^-function or by the integration first by one variable 
e.g. Rl^. Performing the change of variables Rj = kTj/LQ, R'j^ = krl/LQ and 
taking into account that R'^^ > one obtains 



l/l V2 



/"OO ^ ^ 

•^0 i=l fc=l 

= ( A)K+-2-i) / p^2?r'^(^ _ r')e-^", (86) 

I/O "'0 



where r = Ti + . . . + t^^, t' = t[ + . . . + r^2_i, Vt = dn ■ ■ ■ t^^ , and Pr' = 
The integral over Vt' has the form 

/■OO 

9{r) = dT[. ■■T[^_,e{r - r( + . . . + r;_,), (87) 

and may be computed by the Laplace transform over variable r. The result 
is 

The same method applied to the integral over Dr gives 

[VrVT'eir - T')e-'^ = f dTg{T)e~'^ f dn ■ ■■t^,6{t - n - ... - r^ J 



OO ^ui+i^2-2 



{Ui + U2-2)\ ^_(.,+,,_i) 



''^(z/i-l)!(^2-l)! (z/i - l)!(z/2 - 1)! • ^ ^ 

Taking into account Eq. (p9D and that (compare with Eq. (pO|)) 

and 

■k/2 [q a _ 

p\-d<p =ii-=d. (91) 

one obtains 

< dUE)d%{E) >= c:i^Uj^^^^^TYE^, (^2) 
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and 

< c/« >= V c:^c;:r2^ (93) 

But this sum equals 6*2,7-2 (see Eq. (^Sf)) and 

(A'jrf]'^-'^ A 

li^^7^^<C(^)>=T7^=^"' (94) 

"^° <^2n-2 lOTT 

which agrees with the mean part of the bootstrap condition (|22|). 



6 Mean value of Green function products 

Problems similar to the ones described above appear also in the computation 
of mean value of products of the advanced and retarded Green functions of 
the type 

< G"^iE,)GliE2) >, (95) 

where G±{E) = G{E ± it) and G{E) is the Green function of an integrable 
model. A typical example is the semi-classical computation of spectral corre- 



lation functions in the Seba billiard [[lT| , |T^ or, more generally, in integrable 
models with diffraction centers. 

For clarity as in the previous Sections let us consider a rectangular billiard 
with periodic boundary conditions. The Green function for this problem has 
the form 

where e„ = fc^ are eigenvalues of the billiard problem, fc„ is the vector with 
components 27r{m/a, n/b) with integers m, n, and A is the area of the rect- 
angle. This Green function permits also an exact representation as a sum 
over all classical orbits connecting points and x = (x, y) 



G+{x,E) = -'- Yl 4'\k^{x + may + {y + nby), (97) 



m,,n=—OD 



and Hq (x) is the Hankel function, G-{x,E) = G*^_{x,E). 

In diffraction problems one is interested in the Green function when 
X —>■ 0. In this limit the Green function diverges and requires a cut-off 
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at small x (i.e. a regularization) . The divergent part comes only from the 
contribution of the shortest trajectory with m = n = 0. Taking into account 
the asymptotic behavior of the Hankel function at small x one concludes that 
the 'renormalized' Green function (|97|) (when \x\ ^ fi and /i is small ) is a 
sum of 2 terms 

G+{E) = -g4E)+gl^%E), (98) 

where (in 2 dimensions) 



and 



-g^(E) = ^Hk^i) - '- (99) 

9T{E) = -lj:H^'\kL,), (100) 
p 

where Lp denotes the length of a periodic orbit defined by 2 integers m 



and n, Lp = J (may + (nby, the summation is done over all m, n except 
m = n = 0, and parameter /x is a renormalization parameter (a cut-off at 
small distances). In the quantization of the Seba billiard the necessity of 
the renormalization is connected with the fact that the (5-function potential 
is too singular in dimensions > 2 and Eq. ( |100| ) defines one parameter self- 



adjoint extension (defined by /i) of a singular Hamiltonian [|TT|. Note that 
< G± >= g±{E) and < g^'%E) >= 0. 

The semi-classical approximation corresponds to A; — ^ oo and the oscil- 
lating part of the Green function takes the form 

ikLp—3ni/A—eLp/2k 

9T{E) = Y. /, ,, ■ (101) 

ImkLr 



^p 



Let 



gUEi, E,) - ^^^^^ Y. (^^ _ g^ + ,^)„.(^^ _ g^ _ ^^Y ■ (102) 

Assuming that all energy levels have the Poisson distribution (i.e. they are 
independent random variables with mean density d) it follows that 



9mn{El,E2) - -j^^^ J 



dde 



(El - e + ie)'^{E2 - e - ie)" 

'[Cn+m-2- (103) 



27rirf(-l)"-i 
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For rectangular billiard d = A/{16tt) and 

g^E,, E,) = -—Azh^^^c:-:^_^. (104) 



8((Ei-E2)A)™+"-i "+™- 



It is easy to check (see below) that the most divergent part of the product 
(pSf) when E2 -^ Ei equals gmn{Ei, E2) 



< G'l{E{)Gl{E2) >^ g^n{Ei, E^), when E2 -^ E,. (105) 

The purpose of this section is to check that exactly the same answer can be 



obtained by using the semi-classical expression (|101| ) for the Green function 
and applying the saddle point method discussed in the previous Sections. The 
differences between the computation of < d^^{EY^ d^^{Ey^ > performed in 
the last Section and that of < G'^{Ei)G"i{E2) > are (i) the absence of factor 
A/27r and (ii) the change of the definition of e e — *> —i{Ei — E2)/2 (and, of 
course, it is necessary to substitute Ui ^ m and 1/2 ^ n). From Eq. (^) one 
concludes that the saddle point method gives 

om+n-l 4 o 

<G""(^i)G"(^2)> = C^-i 2 r-7^ Z ^ . (^)"^" 

+ v 1^ -V 2j '"+™-2^67r(i(E2-Ei))"+"-i^ A^ 

= *(-l)""'C;i-2^p^3^^^y^4p^T^, (106) 

which agree with the value (|105D obtained by direct calculations. 

Till now we have considered only the dominant term of the mean value of 
Green functions product (|95|) when E2 — > Ei (i.e. we look for contributions 
with the highest negative power of the difference of energies) . But this re- 
striction is not necessary and other terms can also be computed by the same 
method. 

The m-th power of th e Green function (^) may be written as follows 

G'^{E) = --V^ — r (107) 



X 



E 



(E - d -h ie)"^! (E - 62 + ze)™2 . . . (^ _ gp + ze) 



where the first sum is taken over all positive integers m, whose sum equals 
m, 777.1 + . . . + 777p = 777, and the second summation is performed over dif- 
ferent energy eigenvalues e^ such that ci < 62 < . . .. Representing G'^{E) 
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in the similar way and taking into account that the connected part of the 
Green function products should not include mean values of individual Green 
functions one concludes that 

min(m.,n) ^ 

< G']^{E,)Gl{E2) >'= E - (108) 

p=i P- 

sr^ f "^! \ ( ^' \ 

^ 2^ I i i T I I i i T I 9mini9m2n2 ' ' ' 9r, 



mi,n 



. \m1\ni2\ ■ ■ ■ 'mp\ J \n1\n2i • ■ ■ Upl 



im2n2 UmpUpi 



where gmn = 5'mn(-Ei, -E2) is given by Eq. ( p.04|) and the summation is per- 
formed over all positive integers rrii > 1, Ui > 1 such that X^iLi fni = m and 
J2^=i i^i = i^- The factor l/p\ appears from the assuming ordering of variables 
Cj. Because of the symmetry with respect to permutations this factor can 
be removed by assuming the ordering of (say) rrii variables and the above 
expression may be rewritten as follows 

min(m,ra) / I 



< G^iE,)GliE2) >'' = E E 



p=l n^„n. V^l!^2! 
^ 9min\9m2n2 ' ' ' 9 




nipUp 1 

where mi < m2 . . . < mp. 

Each grrnni is proportional to ((£"1 — £"2)^)^""^'""% therefore terms with 
fixed p in these sums will be proportional to {{Ei — £2)^)^"™"" and terms 
with higher values of p correspond to the expansion of < (j™(-Ei)G*" (£2) >'^ 
{El — E2)'^^'^ on positive powers of {Ei — E2)A. 

The above expression represents the value of the connected part of the 
mean value of Green function product calculated from the assumption of the 
Poisson distribution of energy eigenvalues. The total answer includes also 
terms coming from the mean value of the Green function itself (see Eq. (pQ])) 



< G^{Ei)G'l{E2) >=Y.CtGlg\f_ < GX{Ei)G^_{E2) >', (110) 

k,p 

and due to Eq ( |108| ) this expression can also be organized as a series on 
decreasing powers of {Ei — E2)A. 

Let us check that this answer may also be obtained by the saddle point 
method. As was demonstrated in Section ^ the saddle point manifold for 
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smooth terms consists of integer vectors Ni and N', such that 

m n 

EM = E^;- (111) 

The dominant contribution discussed in the previous Sections corresponds to 
the integration over all possible deviations on this manifold. But by fixing 
certain deviations it is possible to find lower dimensional manifolds and the 
integration over them will give corrections to the dominant result. 

Let as above n^ > 1 and mj > 1 be partitions of n and m into positive 
integers: n = ni + . . . + Up and m = rrii + . . . + nip. Each 

such partition gives rise to a possible regrouping of integer vectors 

mi rii 

EN. = EN, 

m2 n2 

i=l j=l 

(112) 

2_^ -i^mi + ...+mp^i+i = 2^ "i+---+"p-l+i' 
i=l i=l 

Such manifold is a part of co-dimension (p — 1) of the full saddle point 
manifold (|1 1 1|) and the summation over all possible deviations on it will 



evidently give S'mmiS'mana ■ ■ -gmpUp- As there exists m!/(mi! . ..nipl) possible 
partitions with fixed rrii and n\/{ni\ . . .Upl) possible partitions with fixed n^ 
the total contribution will be equal exactly to Eq. (|109| ) which demonstrates 



that all correction terms can be computed correctly by the saddle point 
method. 

Let us note that in problems considered there exist 4 parameters with 
dimensions of energy. One is the energy itself or more precisely the center 
of energy window which defines mean values, the second one is the width of 
energy window, the third is the difference of 2 energies, e = Ei — E2, and the 
forth is the mean distance between quantum levels equal 1/d where d is the 
mean density of levels. In the end of this Section we discussed the corrections 
corresponding to powers of the most important dimensionless parameter de. 
All other dimensionalless parameters are assumed to be large and we con- 
sider only the first terms of the expansion on them. The computation of 
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higher order terms of these parameters, coming e.g. from corrections to the 
Gutzwiller trace formulas, are also possible but is beyond the scope of this 
paper. 

7 Correlation functions of classical actions 



It is well known ||13[ , fl^ that there exists a duality between quantum spectral 
statistics and statistics of classical actions. If one is known the other can be 
computed by the Fourier transformation. 

In this Section we demonstrate that the results of the previous Sections 
can be rewritten in the form of correlation functions of classical actions. Let 
us consider for simplicity the 2-point correlation function of periodic orbit 
length for rectangular billiard (calculation of other correlation functions will 



be discussed elsewhere [fH]) 



R^{s) = Y,S{s-l + l')a{l,l'), (113) 

1,1' 

where a{l, I') is a certain weighted function which should dominates by large 
values of / and /'. A typical example is 

a{l,n=6{L-i±^) (114) 

and L is assumed to be a large quantity. The correlation function with this 
weight will be denoted by R2{s, L) 

R,{s, L) = Y.6{s-l + l')5{L - ^). (115) 

i,v ^ 

It is this function which appears naturally in many applications (see below). 

The summation in ( |113| ) is done over pairs of periodic orbit lengths. We 
implicitly assume that all values of / and /' are permitted. There are few 
other possibilities e.g. one can choose /' > /, or only /' > / which lead to 
slight modifications of the formulas below. 

One has 

Ms) = ^ r dTf{T)e'^\ (116) 



27r 
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where 

/(r)=5:e-Ma(/,r). (117) 

1,1' 

This sum sphts into 2 parts 

/(r) = /('^»^)(r) + /("-)(r), (118) 

where 

f''^'^ (r) = Y.ailJ')^p f lail, l)dl, (119) 

1=1' •' 

where we have used the mean density of periodic orbits ( p6D with p = 7c/2A 
and 

/(--)(r) = ^e'^(''-')a(/,/')- (120) 

i^i' 

To compute this function we shall use the saddle point method developed in 
the previous Sections. As above the dominant contribution will come from 
the saddle point manifold ( ^41) and 

/^°''^(^)=EE fdtexp{zTlo{r2-n)+zQt^{---))a{nlo,r2lo), (121) 

where Q = A'^t/{21q). Performing the integration and taking into account 
that ri,r2 —>■ oo but r2 — ri = r is fixed one obtains 



n 



where 



Because 



we finally get 



lo,r V Q\^\ 



D = (iriria(ri/o,Ti/o) = 72 / ldla{l,l). (123) 

(124) 



TT 1 Ap 



Q\r\ ll 2iid^J^¥rl^ 



fir)-/i^^U-J^^c..y, (125, 
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where 

a = ^a(/,/) =p /'/a(/,/)d/. (126) 

I •' 

But the sum in Eq. ( p,25| ) is exactly the semi-classical expansion of oscillating 
part of the trace formula for rectangular billiard (see Eq. ([T^)) with energy 
equals r^ which means that 2-point correlation form-factor for periodic orbit 
lengths is proportional to the semi-classical trace formula 

fir) = ^dir')- (127) 

In particular for rectangular billiard with weighted function ( |114| ) 

4:71 r 
R2is, L) = —L / dTd{T^)e"' (128) 



These results can be considered as a generalization of results of Refs. [|1^ 



[Q . In these papers only the smooth part of this relation has been considered 
/(r) = (i(r^). But relation (|127D contains more information. In particular it 
states that the Fourier transform of classical 2-point correlation function of 
periodic orbit lengths for rectangular billiard should have peaks exactly at 
eigenvalues of quantum problem with Newmann boundary conditions. The 
generalization of this relation to general integrable models is performed in the 
next Section. Numerical calculations for different integrable billiards confirm 
well this prediction |TB . 



The knowledge of the 2-point correlation function of classical actions per- 
mits to compute easily any quantities depending on a product of 2 densities. 
For example, 

A ik{l-V)-e{l+l')/2k ja .oo , ^r 

dl{E) = 2i—Yy- = = ^^ dsd^'-^'^'^R^is.L)^ 

= 2S^(^)i e--"^dL = -d(E), (129) 

which coincides with the bootstrap condition (pl|). 

It is also instructive to compute the 2-point correlation function of energy 
levels 

R,ie) = < diE + ^)diE - ^) > 

A ik{l-l')+e{l+V)/ik 

= < d\E) > +{-f <Y. ^^^ +c.c.>. (130) 
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The last sum equals 

(131) 

and the 2-pomt correlation function for energy levels of rectangular billiard 
(and for general integrable systems as well) has the following form 

i?2(e) =< d\E) > +6{e) < d{E) > . (132) 

Here the brackets denote a local average over a certain energy window. When 
this window is very large < d{E) >= d but in general this smoothing defines 
a new local scale and < d{E) > can deviates from d. This result is exactly as 
was expected. The energy levels of typical integrable systems (in particular 
for rectangular billiard) are distributed as independent random variables but 
at the local scale (or equivalently only after unfolding). The correlation 
function of classical lengths ( |127| ) does the required transformation from 
global scale to the local one. 

8 General integrable systems 

In the previous Sections for clarity we discussed the case of rectangular bil- 
liard but the method used is not restricted only to this example and can be 
generalized for general integrable systems as well. 

The starting point of the derivation of the trace formula for integrable sys- 
tems is a formal expression of a Hamiltonian as a function of action variables 

E = H{h,h). (133) 

Semi-classical quantization conditions consist in fixing the values of these 
action variables |]T|-[^ 

Ii = n{ni + -fii), (134) 

where rii are integers and /i, are certain phases connected with the type of 
classical motion (the Maslov indices). Therefore for 2-dimensional integrable 
systems (multi-dimensional case will be discussed elsewhere [0) the quan- 
tum energies are 

1 1 

-Enma = H{ni + -Hi, n2 + -/i2), (135) 
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and the quantum density of states is 

d{E)= ^ 6{E-E,,,^,). (136) 

ni,n2 

Using the Poisson summation formula one obtains 

d{E)= J2 /e'"'(^^"^+^^"^)<5(E-E„,„Jrfnirfn2. (137) 

Ni,N2 

Separating the term with Ni = N2 = and setting rii = li/h — fii/4: leads to 

d{E) = d{E) + d'"%E), (138) 

where 

d{E) = ^J5{E- H{h, h))dhdh = 7^^ / 5{E - H{p, q))dpdq, (139) 
and 

d°-(E) = \ J2 e-'(^i^i+^2^2)/2 f e^^'^''^'^+''^''^/^5{E - H{h,l2))dlidl2. 

(140) 
It is convenient instead of variables / to use variables e and t defined by 

where the unit vectors v and r have the following components 

1 1 



-LOi, Ti 



J2e^,Uj, (142) 



and 

0.. = -^j-- (143) 

are frequencies of the classical motion, and Cij is 2 x 2 antisymmetric tensor 
with components en = 622 = and 621 = —612 = 1. 

Vector u is the unit vector perpendicular to the line of fixed energy, vector 
r is the unit vector tangent to this line, and z7^ = f^ = 1, z7r = 
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It is easy to check that dlidl2 = dedt and 



J5{E-H{e))de= } (144) 



-'2 

The integral over the surface of constant energy 

^2.i{N,i,+mi,)in^^ (145) 

can be taken by saddle point method. The saddle points are points where 

EA^^f = 0. (146) 

This condition means that in saddle points the classical frequencies should 

be commensurable 

ujx _ Ni 

UJ2 N2 

Straightforward calculation gives 



(147) 



where the curvature K is 

As 

d^Nih + N^h) ^^r^n K 2 



dt^ fr[ dt 

the oscillating part of the level density takes the form 






d"'^\E) = ^ J2 Pexp{l-S,i{E)-t^{Nifii + N2fi2)-tjSgnK)+c.c., 
'^ Ni,N2>o a z 4 

(151) 
where the pre-factor P is 

P= I ^ ^^- (152) 

l{NiuJi + N2UJ2)\K\ 

30 



The classical action of resonant periodic tori (= periodic orbit) 

Sa{E) = 2Ti{NJ^ + N2l2), (153) 

is defined by special values of action variables Jj = Ii{Ni,N2,E) for which 
the following two equations are valid 



and 



H{h,h) = E 



^ mh,h) ^^mh,h) 



dh dh 

Introducing the classical period of motion 

^ ' dE cui' 

one can rewrite the pre-factor in the following form 



(154) 
(155) 

(156) 



where 



K 



Nl 






P 



+ Ni 



TjE) 



d^H 



2N1N2' 



d^H 



(157) 



:i58) 



Eq. ( |151| ) is the semi-classical trace formula for general integrable systems 
and is a generalization of the trace formula for rectangular billiard given by 
Eq. (|T^. The important difference between these two formulas is that the 



periodic torus action is given by a simple formula ([lID for the latter but only 
by implicit formulas ( |15iJ| )-(l55) for the former. 

The first step to apply the saddle point method similar to the one dis- 
cussed in previous Sections to the general trace formula (|151|) is to find the 
expansion of the classical action Sci{Ni + 6Ni, N2 + SN2, E) into series of 6Ni. 

By differentiating Eqs. (|153|) -( pr55D one obtains that this expansion (which 
is a generalization of Eq. (ESI)) has the following form 



S,i{N, + 6N^,N2 + 6N2,E) = S,l{N^,N2,E) + 27^{h6N, + I26N2) 

+ Q{N26N, - N,6N2)\ (159) 
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where 

and all other notations are the same as above. 

Taking into account that /i,/2 do not depend on the common factor of 
integers A^i , N2 and assuming that values of action variables /j for different 
primitive periodic tori are non-commensurable one concludes that the only 
possibility of existence of non-trivial saddle points (i.e. the cancellation of 
linear terms from different periodic tori in the exponent) is exactly the same 
conditions as in previous Sections. For example, for two periodic tori the 
saddle point manifold has the form exactly (up to notations) the same as in 



Eq. (H) 



Arf ) = r(2)ni + 5Ni, N^ = r^^^n2 + SN2, (161) 

where superscript denotes values for different periodic tori. As before, we 
assume that integers rii and ^2 have no common factor. 

All further calculations can be performed exactly as it has been done 
for rectangular billiard in Section 4. The only difference with the latter 
case is different expressions of the pre-factor in the trace formula ( |157| ), the 
coefficient of the quadratic form (|16CI|) , and the period of the motion. But it 
is easy to check that 

which agree with Eq. (pO|) and all formulas derived for the oscillating parts of 
products of periodic orbit contributions for rectangular billiard remain valid 
for generic integrable systems as well. Note that from Eq. (|45|) it follows that 
the total Maslov phase will have the correct value. 

To apply the saddle point method to smooth terms one need to know 
the mean density of periodic orbits for general integrable systems. Though 
it is possible to do it by applying the Hannay-Ozorio de Almeida method to 
integrable systems 0] it is instructive to compute it directly from Eqs. ( |153|) - 



Let us define the density of periodic orbits with fixed period ( [156| ) as 

POD 

p{r) = / 6iT{Ni, N2, E) - r)dNidN2. (163) 
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Introducing polar coordinates Ni = R cos (f), N2 = R sin one gets 
,27ri?cos< 



^ r 27r^Rcos0 _ ^^^^^^^ ^ r r ^ ^^^^^ ^^g^^ 



This integral can be transformed into an integral over classical actions Jj by 
noting that they are functions of and E defined by Eqs. ( p.54| ) and ( |155| ). 
Straightforward calculations give 

det ( ^^'1^^ ^^'1^^ \ - 1 (165) 



where K is the curvature defined in Eq. ( |149D . As 002/ uJi = tan 



n2 rr 3°^ H cP'H 

K = cos^ ^~qTT + ^^^^ ^~dP' ~ ^ ^°^ *^ ^^^ ^'dTdT' ^^^^^ 

Therefore 

d(j) = K5{E - H{h, I2))dhdl2, (167) 

and in particular 

p{r) = -^ J{ujI + ujI)K6{E - H{Iu I2))dhdl2. (168) 

The computation of smooth terms can be performed exactly as above but 
with different values of pre- factor, P and the coefficient, Q, of the quadratic 
form in the exponent. We have 

and 

Q = ^. (170) 

COS0A 

Because the reduced period 

T = 2vr^ (171) 
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the combination of these terms in Eq. ( |162| ) has the same value as before and 
all integrals except the integral over rest exactly as in Section ^. The last 
integral takes the form 



J Pj^dcp = J^ = Js{E- H{h, h))dhdh = d, (172) 



as it should be to reproduce the smooth part of bootstrap conditions (]22|). 

Though one can compute any correlation functions of classical actions 
with arbitrary definition it is convenient instead of Eq. ( |113| ) to define the 



2-point correlation function of differences of actions for general integrable 
systems as follows 

Ms) =Y.^is-Sp + Sp,)P,Pp,e-'^^^^"^-^'^a{Tj,, T^,), (173) 

p,p' 

where Tp are periodic orbit periods, /ip are the Maslov indices for the problem 

considered, Pp are pre-factors in the trace formula given by Eq. ( |157D (without 

powers of the Planck constant), and a{Tp, Tpi) is a certain weighted function. 

Separating the diagonal and oscillatory contributions as in Section ^ 

R2{s) = ^ f dTe'^'{f'~''^^\T) + j<"''\t)), (174) 

and using Eqs. (|169|) -( [l72l) one gets 



where 



f^'^'Hr) =j:P>iTp,Tp) = ^ J dTa{T,T) J ^ = ad, (175) 

a = — f dTa{T,T). (176) 

The computation of oscillatory terms can be performed exactly as in Section 
1^. Taking into account the relation (|162|) it is easy to check that 

j:[osc)f\ ^ _^ y- p ^iTSp-i7rfj.p/2~i7rSgnK/i ^ Q-j^^N 

V^ p 

Therefore 

fir) = n''diE,h)\n=yr, (178) 

where d{E, K) is the level density of the quantum problem considered which 
generalizes the results in Refs. [0 , [|14| about duality in integrable systems. 
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9 Conclusion 

The main goal of this paper was to develop a new method which permits to 
calculate sums of product of periodic orbit contributions for generic integrable 
systems beyond the diagonal approximation. In the previous Sections it was 
demonstrated how different quantities like bootstrap conditions (^), mean 
value of Green function products (p3), 2-point correlation function of classical 



actions (|113|) , etc., can be computed for integrable systems directly from 



semi-classical trace formulas using this method. More applications of this 



method will be discussed elsewhere 16 



The main ingredient of the method is the existence of hidden saddle points 
in the summation over periodic orbits. The saddle point configurations con- 
sist of periodic orbits whose integer vectors are mutually parallel. The saddle 
point manifold includes orbits which (i) are almost mutually parallel and (ii) 
are chosen in such a way that the total action is quadratic on deviations from 
the saddle point configuration. The integration over these quadratic forms 
gives the dominant contribution to the summation over periodic orbits. By 
considering low-dimensional sub-manifolds it is also possible to compute cor- 
rections to this term. 

In the calculations it was assumed that for generic integrable systems 
classical actions of (almost) all resonant tori are non-commensurable but the 
method can be generalized for systems with geometrical symmetries and (or) 



accidental degeneracy. In Ref. [|T^ the method has been applied for integrable 
quantum maps. 
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